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Exact Chiral Symmetry on the Lattice: QCD Applications* 
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I review recent progress and results in lattice QCD obtained using fermions with exact chiral symmetry. 


1. Introduction 

In the last few years it has been understood 
that exact chiral symmetry can be realized at fi¬ 
nite lattice spacing. A breakthrough started with 
the observation that massive Dirac fermions in 
4+1 dimensions reduce to chiral fermions in 4 
dimensions under certain conditions [:l],||||]. The 
resulting effective operator D of light boundary 
fields |4|i| has the correct continuum limit, no 
doublers, and is local |J. Most remarkably, it 
satisfies the Ginsparg-Wilson (GW) relation (Q 

75D + £>75 = 0D75D , (1) 

which guarantees an exact chiral symmetry [|| 

= 75(1 - aD)q , 5q = q^ 5 (2) 

of the fermion action at non-zero lattice spacing. 

Within the perfect-action approach ||, a 
fermion operator that satisfies the GW relation 
can be defined but no explicit construction 
has been found so far (see |tl| for a recent review). 

In QCD a chiral symmetric regularization en¬ 
tails many theoretical advantages: in the infrared 
it allows one to simulate massless quarks and it 
provides a natural definition of the topological 
charge; in the ultraviolet it simplifies the subtrac¬ 
tion of divergences in composite operators. 

Last year an important step forward was the 
demonstration that Neuberger’s fermions can be 
used for large-scale QCD computations, at least 
in the quenched approximation. As a result, very 
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light quarks can be handled and legendary prob¬ 
lems such as the A/ = 1/2 rule in K —* nn decays 
are greatly simplified and can be attacked. 

In this talk some of the phenomenological ap¬ 
plications of Ginsparg-Wilson fermions are re¬ 
viewed. Most of the numerical studies reported 
are exploratory. Please refer to the parallel ses¬ 
sion contributions for interesting theoretical de¬ 
velopments which are not covered in the follow¬ 
ing. 


2. Domain-wall-overlap fermions 

The five-dimensional domain-wall Dirac opera¬ 
tor can be defined as 


V = 


1 

2 L 


j 5 (d* + d s ) - a s d*d s 


X 


( 3 ) 


where s labels the sites and a s the lattice spacing 
in the fifth dimension. The operators d* and d s 
are the forward and backward derivatives, 


X = D\y — — 
a 


a 

a = — 
P 


( 4 ) 


with 0 < p < 2. The massless four-dimensional 
Wilson operator is 


Dw = 


1 

2 L 


7m(Vm + V;)-oV*V, 


( 5 ) 


where and V* are the gauge-covariant forward 
and backward derivatives 


V^g(x) = ^ U^(x)q(x + afi) - q(x) 

V*g(x) = - q(x) - Ufa - afi)q(x - ajl) 


(6) 
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and a is the lattice spacing in the physical four 
dimensions. The operator is supplemented with 
the boundary conditions 


P+q( 0, x) = P-q(a s N s + a s , x) = 0 


(7) 


where P± = ^(1 ± 75 ) and N s is the extension 
in the fifth dimension. From a four-dimensional 
point of view, the system corresponds to QCD 
with many flavours mixed in a peculiar way. 

By integrating out the heavy flavours, a four¬ 
dimensional effective action of light boundary 
fields is left 


N, 


aD Ns = 1 

where 
X = 


■75 


ClsX 


(1 + H) N ‘ - (1 -H) 

(1 + H) N * + (1 - H) N ° 


H = 75 X . 


2 + a s X 

For N s —> 00 , a massless effective action 


ciDdw = 1 + X- 


1 


Vx^x J 

is derived and the Neuberger operator || 

clDN = (1 ~\~ X 


( 8 ) 


(9) 


( 10 ) 


( 11 ) 


VxTx, 

is obtained if the limit a s —> 0 is also taken. Most 
remarkably, the effective massless Dirac operators 
in Eqs. © and ( |lT| ) satisfy the GW relation Jl3|. 


3. Exact chiral symmetry at finite a 


For a given operator D that satisfies the 
GW relation, the QCD fermion action with Nf 
flavours can be written as 


Sf 


= \(D + P-M ] P- + P + MP+)i/j 1 Or) 


where 

-P± = ^(l±75), 75 = 75(1 ~aD), (12) 

M = diag(mi,..., m Nf ), ip = (q 1 , ■ ■ ■, q N f ) and 
ip is defined analogously. It is invariant under the 
U(-/V/)lxU(1V/) r global transformations 

Ip-L VlV’L i’L -> 

4>R -> VrV'R ^R -> "OrDr , (13) 


where Vl,r eU(7V/) L , R and 


V’R.l = P±i> V’R,l = ipP y , (14) 

if also M —> No additive quark mass 

renormalization is required. The action is O(a)- 
improved, since no chiral invariant operators of 
dimension d = 5 can be constructed. 

The global chiral anomaly is recovered a la Fu¬ 
jikawa The fermion integration measure is 

not invariant under U(1)a transformations, and 
the topological charge density from the corre¬ 
sponding Jacobian 


a 


a Q{x) = — Tr 7 5 D(x,x) 
2a 

satisfies [51151 


(15) 


n_ — n+ = index(D) = a 4 Q(x ) , (16) 

X 


with n± the number of right and left zero modes 
of the fermion operator. 

Bilinear fermion operators with proper chiral 
transformations 

0^(x) = q a {x)TqP{x ), qP =(l- ^D^jqP (17) 

are 0 (a)-improved, but they do not transform in 
a simple way under CP |l6| , |17| . However in cor¬ 
relation functions of local operators at non-zero 
physical distance, it holds 

= n _L„ ( 18 ) 

(1 2 17l l 3 > 

and a simple CP transformation [x = (xq , —x)] 

Of(x)^ 1 -~P^O 0 r a (x) (19) 

1 2 m 0 

is recovered ||l 6 | . The generalization to four- 
fermion operators is straightforward. 

Non-singlet local rotations lead to exact vector 
and axial Ward identities (WI), and the very same 
definition of the bare quark mass appears in the 
axial and vector WIs and in the quark propagator. 
The conserved currents V“ and A 1 !' can be con- 

P P 

structed by extending the gauge group SU ( N c ) —► 
SU(N c ) x U (1) M. By performing a local U (1) 
flavor rotation 


U^x) - U^\x) = e ia - 


{x) U^x) , 


( 20 ) 
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the kernel 


K fl = -i 


.SD(U^ a) ) 


Sa^(x) 


a=0 


( 21 ) 


can be used to define the following conserved cur¬ 
rents [fist] 


V“ = g(P-V + +W_)r g (22) 

a; = g -(?_v + -p + v_)r 9 ( 23 ) 


where T a is a generator of the non-singlet trans¬ 
formations. It is interesting to note that (con¬ 
served) current-current correlators and their gen¬ 
eralizations require the propagator from any 
point to any point (h]] . 

In the chiral limit, the local anomalous flavour- 
singlet WIs read 

(d;A^(x)d) = 2Nf(Q(x)(D) + {$%&> , (24) 


where 8 a O is the local variation of any finite 
(multi)local operator O, and A^x) is the sin¬ 
glet axial current. By assuming the absence of a 
(7(4(1) massless Goldstone boson, the correspond¬ 
ing integrated WIs read 


0 = 2 Nfa A ^{Q{x)6) + (8 A 6) . (25) 

X 

Since the second term in the r.h.s. of Eq. (|25| ) is fi¬ 
nite, it follows that a 4 Y^ x Q{ x )is also finite, as it 
has finite insertions with any string of renormal¬ 
ized fundamental fields. Therefore Q(x) can only 
mix with operators of dimension < 4 and vanish¬ 
ing integral, hence only with d*A^(x). No power- 
divergent subtractions with lower dimensional op¬ 
erators (such as the pseudo-scalar quark density) 
have to be performed |il|. This is a very dis¬ 
tinctive feature of GW fermions. Calling Z the 
mixing coefficient, one can define finite operators 
Q and .A” by writing 

Q{x) = Q{x) - ^-d*A 0 ^(x) (26) 

Al(x) = (l-Z)A°^x), (27) 


and the renormalized singlet axial WIs are 
(d*Al(x)6) = 2N f (Q(x)d) + (S\6) . (28) 


Renormalization constants of several compos¬ 
ite operators have been studied at one loop in 
perturbation theory for domain-wall fermions 


21 

BO 


22.23. 
31 


25| and overlap fermions 


20 


2 !) 


. Non-perturbative determinations have 
been obtained for the non-singlet local vector and 
axial currents MEOEL. for the scalar and pseu- 

36], and for some four- 


doscalar densities 4B3 34 


fermion operators |57 38 


4. Meson spectroscopy 


In the past year several collaborations have 
computed the meson spectrum and light quark 
masses using GW fermions, either with the over¬ 
lap formulatio n f^ , j35|]3G| , |3i| or with the per¬ 
fect actions [ |lo|j4lf . Simulated lattices have lin¬ 
ear extensions L = 1-3 fm and lattice spacings 
a = 0.08-0.2 fm. Results for the pion mass 
squared in units of I/j-q (ro = 0.5 fm) as a func¬ 
tion of the quark mass normalized at the refer¬ 
ence point AIp = 2 M\ ( Mk = 495 MeV), are 
shown in the first plot of Fig. [j]. Data in the range 
500 < Mp < 800 MeV show a linear behaviour 
with a vanishing intercept within the statistical 
errors. The linearity manifests itself as a wide 
plateau in the second plot of Fig. [I] and it is in 
very good agreement with what was previously 
observed with Wilson-type fermions in the same 
range of masses (see for example HSU)- For 
degenerate quarks, quenched chiral perturbation 
theory at the next-to-leading order (NLO) pre¬ 
dicts [fdsl.fdfj 


M 2 


E 

2 to F 2 . 

aM' 2 


1-6 


3(4t rF 1 ) 5 


(—(f))- 

(l + 2 l0g(^)) + (2„- Os ) 


(29) 


M 2 


(4ttF) 2 J 


where E, F, mo, a are the leading-order (LO) 
couplings of the quenched QCD chiral Lagrangian 
|^,Q, « 5 ,a 8 are some of the NLO ones, p x = 
4t tF, M 2 = 2Em/F 2 and 8 = m$/3(47rF) 2 . A 
comparison of Eq. (^9|) with the data in Fig. |l] 
indicates small corrections due to quenched chiral 
logs and/or higher order terms in this range. 

Since GW fermions do not suffer from excep¬ 
tional configurations, lighter pion masses can be 
simulated, provided the physical volume and the 
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• GHSW 

• GHR 

• BGR Coll. 


m/m. 


Figure 1. On the left, (ro-Mp) 2 vs the quark mass normalized at the reference point Mp = 2 M\\ green 
circles ©I, black circles Q, blue squares ©, green diamonds 63], red circles [ph yellow diamonds 
f39f . On the right, (roMp) 2 /(m/m re i) vs (m/m re f): green circles |4j| , black circles [pH , red circles pdf . 


cut-off are large enough. First exploratory stud¬ 
ies in the region 200 < Mp < 400 MeV have been 
reported at this conference |3| US 



Figure 2. y vs x with FP action from Ref. fill. 


The presence of quenched chiral logs can be 
tested in the double ratio of meson masses with 
non-degenerate quarks 147 44| ] 

4?nim2 


V = 


Mp 12 Mp 12 


(mi 


m 2) 2 Afpn Mp 22 


(30) 


At NLO 


y = 1 + Sx + + 0(m\,rnl) , (31) 

where 


x 


w 


mi + m 2 


lo g p 

\m-\ 

( MiT7i2 , (m 2 

- log I - 

\mo — m 1 \m i 


m 1 — m 2 
2mim 2 


mi ) 
)- 


m 1 



(32) 

(33) 


and p x appears explicitly in Eq. (|3^) at higher 
orders only. 

Using the fixed-point (FP) and chirally im¬ 
proved (Cl) actions, the Bern-Graz-Regensburg 
(BGR) collaboration computed meson masses 
with non-degenerate quarks In Fig. |2j 

data obtained with the FP operator on a lattice 
with a ~ 0.15 fm and L ~ 2.4 fm are shown. 
By neglecting 2aEw/3F 2 (47rE) 2 and higher or¬ 
der terms in Eq. (©), they obtain <5 = 0.17(2) 
and 5 = 0.18(2) for the FP and Cl actions respec¬ 
tively ©]. More studies are needed to properly 
assess the systematics due to finite volume effects 
and/or finite lattice spacing, and to remove the 
uncertainties due to the leftover explicit symme¬ 
try breaking of these actions. An extensive com¬ 
parison of these results with those obtained in the 
past with standard actions can be found in |5(J . 
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5. The chiral condensate 

The chiral condensate E has been computed 
with Wilson-type fermions by fitting Mp/2m 
with the LO term on the r.h.s. of Eq. ( |29| ) 
pTlp^ . The same analysis has been repeated 
with overlap fermions in Refs. |34p(| and the •§, 
results are reported in Fig. ||| as well. Within w 
the statistical errors, the results are compati¬ 
ble with the continuum extrapolated value ob- ^ 
tained with non-perturbatively improved Wilson c« 
fermions p3| . This represents a further indica¬ 
tion that 0(a 2 ) effects are moderate with overlap 
fermions. More studies with lighter quark masses 
and larger volumes are needed to properly assess 
the systematics of these encouraging results. 
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Figure 3. Chiral condensate vs a: black squares 
pill , black diamonds ||52|, black circles |h|], red 
circles Q, red squares H- 


Properties of QCD Green’s functions in a fi¬ 
nite box of linear extension L and with very light 
quark masses can be studied within chiral pertur¬ 
bation theory J54|. If 2mSL 4 ~ 1, L> 1/(47 tF) 
and p 2 ~ 1/E 2 , the correlation functions can be 
expanded in powers of a parameter e with 


•\/ 2 £ m 

A 2 

x 



(34) 


Figure 4. Quark mass dependence of the 
scalar condensate for three volumes: 8 4 (circles), 
10 4 (squares), and 12 4 (triangles), from Ref. p3| . 


and A x being the cut-off of the effective theory 
|54[ |. At leading order the partition function is 
given by 


Z(m, 9) 

rl 

ex P 77 


dUod £ x 


(35) 


JSV(N f ) 

d 4 xTr (d^d^) + zReTr ( e i9 / N fU 0 ) 


where the pion field is factorized as U{x) = 
Uq exp (i\/2£{x)/F) and z = mSF. The integral 
over the global mode 


/ dUo exp 

Jsu(N f ) 


zReTr {e i9 / N fU 0 ) 


(36) 


needs to be done exactly, while a reordered chiral 
perturbation theory applies to the non-zero in¬ 
tegration mode £(x). Partition functions Z v (m) 
and correlations can be defined in sectors of fixed 
topology by Fourier transforming in 8 J55|. By 
comparing the chiral perturbation theory expec¬ 
tations for the correlation functions with the lat¬ 
tice data, the basic assumption of spontaneous 
symmetry breaking can in principle be verified 
and the low energy constants (LEC) of the chiral 
Lagrangian extracted. Properties of QCD in the 
infinite volume limit can then be recovered. 
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The e-expansion has been extended to 
quenched QCD in Refs. US At the leading 
order, the chiral condensate in a sector of fixed 
topology v reads 

= z(^I v (z)K v (z)+I I/+1 (z)K u _ 1 (zfj+^(37) 

where I v {z) and K v {z) are modified Bessel func¬ 
tions [|3 ■ One-loop corrections give |58| , [59| 

Si- lo °P(z) = —T, v {z') (38) 

z 

where z' = mE e ffV, 




1 - 


a 




3(4t tF} 


r(A + *(zl)) 


3(4t tFL) 


T01 


(39) 


1/To is the renormalization scale and /3, are 
two universal “shape coefficients” J5f|. It is in¬ 
teresting to note that the expansion parameter 
e ~ l/(4irLF) in Eq. ( |38| ) is comparable to 
Mp/(4irF) 2 in Eq. (|2|) for light pions (Mp ~ 
1/L), while the prefactors turn out to be differ¬ 
ent. 

Quenched QCD in the e-regime has been ex¬ 
plored on the lattice in Refs. 60,53||. In Fi g. |~i| 
an example of the results obtained in Ref. |53| 
with overlap fermions for the vacuum expecta¬ 
tion value of the scalar density in the topologi¬ 
cal sector v = ±1 at different volumes and var¬ 
ious masses is shown. Even with poor statis¬ 
tics and small volumes, a signal compatible with 
the expectations of chiral perturbation theory has 
been reported. Analogous analyses have been at¬ 
tempted in Ref. |n],[io). More studies at larger 
volumes and higher statistics are needed to con¬ 
firm the indications of these exploratory investi¬ 
gations. 


6. Topological susceptibility 

In the chiral limit, the Fourier transform of the 
singlet axial WI in Eq. (^8|) for O = Q reads 

Xt{p) = a 4 ^e-^(Q(a;)Q(0)) + CT(p) (40) 

X 

= ^ E (x)Q(O)) + CT(p) 

f x 


where the same contact term CT(p) has been 
added to both sides of Eq. (|t]) to make them 
separately finite. CT(p) is a fourth-degree poly¬ 
nomial, which can be chosen to vanish at p = 0 
since in this case the second line of Eq. ( |40| ) is 
certainly finite in the absence of zero-mass par¬ 
ticles in the singlet channel. As a consequence, 
the topological susceptibility vanishes in the chi¬ 
ral limit, i.e. Xt(0) = 0 [|l9f . 

For Mp —* 0, the leading chiral behaviour 

rp2 71^2 

X t (0) = ^ + O(M 4 ) (41) 

is recovered by comparing the integrated singlet 
axial WIs with O = Q and O = P° j6/| . 
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Figure 5. Summary for the quenched topological 
susceptibility from [|l0|: filled boxes Q, empty 
triangles Q , filled circles |G7j, empty circles Q , 
filled diamonds Q, filled triangles p9| . 


In the chiral limit, Xt(p) satisfies a three-times- 
subtracted dispersion relation 


Xtip) = h + b 2 p 2 + b 3 (p 2 ) 2 
?2 

(P 2 ) 3 




p(t ) 


p 


M 2 


(t + p 2 )t 3 


(42) 


dt 
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where the contribution of the rf meson has been 
separated since it is expected to dominate the dis¬ 
persive integral 163 . For p 2 —> 0, the “sum rule” 
Xt(0) = 0 implies 


F 2 ,M 2 , 

, _ 77 77 

1 21V f 


(43) 


where R 2 , = F 2 ,M^,/2Nf. Under the “smooth¬ 
quenching hypothesis”, the Witten-Veneziano 
formula is obtained p4| , |65f 


F 2 M 2 

Wf 


-=o 


(Q{x)Q{ 0))ym 

X 

((n_ — n+) 2 ) 

V 


(44) 


where the quantum average has to be done in the 
pure Yang-Mills (YM) theory JuJ. 

By using Wilson and staggered fermions it was 
argued that [ [toIJtiI 

Xt( 0) = ^r o (^) 2 a 4 ^(P 0 (a:)P 0 (0))^ en (45) 

J X 


where only Zweig-violating (ZV) diagrams are in¬ 
cluded. With GW fermions, this formula is the 
algebraic equivalent of Eq. (|44|) |b^| . 

The topological susceptibility defined by us¬ 
ing Neuberger’s operator has been computed for 
several lattice spacings and volumes for a pure 
SU(3) YM theory |6f|. Analogous computations 
have been performed in the past year with over¬ 
lap fermions [|69| , FP and Cl actions p0| , |66|| . A 
summary of the results obtained can be found in 
|40| and is reported in Fig. ||. The central value 
of Xt obtained with GW fermions is quite stable 
as a function of the lattice spacing and is also 
compatible with the value obtained by other ap¬ 
proaches j72],[ 73 ]]. More work is needed at larger 
volumes and smaller lattice spacings before these 
encouraging indications are fully confirmed and 
the magnitude of the systematic error is properly 
assessed. 

Exploratory computations of \t f° r Ac > 3 
have been performed in the past year [[7~i| . |77| . [7b| . 
The results are compatible with a smooth large- 
N c limit and a non-zero \t i n the limit N c —> 00 . 


7. K-> 7T7r decays 


Non-leptonic I\ 
parametrized as 


7T7T amplitudes can be 


T[K+ -> 7T+7T 0 ] 
T[K° -> TT+TT"] 
T[K° -*■ 7r°7r°] 


(46) 

\[\ A 0 e iS ° - 2\[\A 2 e lS * 


where <5/ and Aj are the 7T7 t phase shifts and the 
isospin amplitudes for I = 0, 2. Direct and indi¬ 
rect CP violation are parametrized by 


1 J$ ReA 2 /Im A 2 InrA 0 \ 
\[2 11 e .1 () y He/1 2 He. do J 


(47) 


T[K l -4 ( 7 T 7 r) 0 ] 
T[K S -> ( 7 T 7 r) 0 ] 


(48) 


respectively. Experimental results reveal 4> = 
7r/2 + 5 2 — <5o ~ 7 t/ 4, a A I = 1/2 selection rule 
|A 0 /A 2 | ~ 22.2 and the presence of direct and 
indirect CP violation in nature: 


Re(e7e) = (16.6 ± 1.6) x 10“ 4 0 

| £ | = (2.282 ± 0.017) x 10“ 3 (7| . (49) 

Phenomenological analyses of the unitarity trian¬ 
gle indicate that the Standard Model picture of 
indirect CP violation in the kaon system is consis¬ 
tent with that of B decays and oscillations |79|,|o| . 

The A I = 1/2 rule and the value of e'/e can 
be explained within the Standard Model only if 
the strong interactions crucially affect these non- 
leptonic weak transitions (see pT| . p2p^ ] for recent 
reviews). In this case a more complicated blend of 
ultraviolet and infrared effects prevented reliable 
determinations of the relevant matrix elements. 
Power-divergent subtractions can be needed to 
construct the renormalized operators that enter 
the effective Hamiltonian [^4|j85| , |8^ , ^l| . In the 
infrared, the continuation of the theory to Eu¬ 
clidean space-time and the use of finite volumes in 
numerical simulations generate a non-simple re¬ 
lation between the physical amplitudes and those 
computed on the lattice 
















7.1. The A I =1/2 rule 

By using the operator product expansion 
(OPE), the CP-conserving AS = 1 effective 
Hamiltonian above the charm threshold is given 

by 


TT AS—1 


Gp 

V 2 


C + ( M )0 + (M)+C'_(/x)a_( M ) 


where the Wilson coefficients C± (/x) are known at 
the NLO |j(|[n]] and the bare operators are 


o± 


(s a 7ll P-U b )(u\P-d a ) 


(50) 


± (sj fi P-u)(u'y fl P-d) 


(u -> c) . 


The contributions that arise when the top quark 
is integrated out are heavily suppressed by CKM 
factors and can be neglected. 0± belong to dif¬ 
ferent chiral multiplets and are CPS-even. In 
correlation functions at non-zero physical dis¬ 
tance, 0± cannot mix between themselves or with 
other four-fermion operators, but only with the 
dimension-six operator |fsT[fi7)| 


Qm = (m£ - m 2 c ) m d (sP+d) + m s (sP-d) . (51) 


The renormalized operators are 


0 ± (ji) = Z ± (jj) 


0± + 6™ Q n 


(52) 


The so-called QCD-penguin operators are 

23,5 = (s7 m P_J) Y ( 55 ) 

q—u,d,s 

Q 4.6 = {s a lti p_d b ) Y (9 b 7(56) 

q—u,d,s 

(see Refs. |8l],|2| for definitions of the other op¬ 
erators). At non-zero physical distance, mixing 
with two lower-dimensional operators 

Q P = m d {sP+d) + m s (sP-d) (57) 

Q c r = mdisF^a^P+d) + msisF^a^P-d) 


can occur and power-divergent subtractions are 
needed even with Ginsparg-Wilson fermions. 

With Wilson fermions, only CPS and flavour 
symmetry can be used to determine the renor¬ 
malization pattern of 0 . Even with an active 
charm, a quadratic divergent contribution needs 
to be subtracted in the parity-conserving sector 


= Z ± M 


ol c 




b T +Q T 


b% 

-4-Qi 


(58) 


where 


where Z±(p) are logarithmic-divergent renormal¬ 
ization constants and are suppressed by a 
factor a s . No power-divergent subtractions are 
needed to renormalize 0± when fermions with an 
exact chiral symmetry are used |nj . 

For m s ^m d , 


Qm 


m c )d* 


- m d + m s 
-m s - m,d 



m d - m s Asd 
m s + m d M ■ 


(53) 


and it does not contribute to matrix elements 
which preserve four-momentum [ ^~i||85|i . 

If the charm is integrated out not only poten¬ 
tially large contributions of 0(pi 2 /m 2 c ) are ne¬ 
glected, but ultraviolet power divergences can 
arise in the renormalization pattern of the rel¬ 
evant four-fermion operators. In this case the 
AS = 1 effective Hamiltonian can be written as 
r 10 

H ^ S=1 = (54) 

V 7—1 


Q s = (m„ - m c )sd (59) 

Qt = K — m c )s a^F^d . (60) 

and are four fermion operators with wrong 
chirality. In this case only the flavour part of 
the GIM mechanism survives due to the explicit 
breaking of chiral symmetry J85| . 

The LECs of the CP-conserving AS = 1 elec- 
troweak chiral Lagrangian can be extracted from 
three-point correlation functions, thus avoiding 
the infrared problem that affects the direct com¬ 
putation of the K —> 7T7T matrix elements on the 
lattice §. No large cancellations among leading 
order terms are expected in the ratio Aq /A 2 |; an 
enhancement should therefore be visible already 
at this order. As a result, a combined use of 
fermions with an exact chiral symmetry and chi¬ 
ral perturbation theory can be the starting point 
to attack the A I = 1/2 rule. 

In the past years the RBC and the CP-PACS 
collaborations have studied the A/ = 1/2 rule 
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with domain-wall fermions with a finite fifth di¬ 
mension |37||j2|]. They have computed the K —> tt 
and K —> 0 matrix elements for the operators of 
the AS = 1 effective Hamiltonian in Eq. (^d|) and 
used LO chiral perturbation theory to recover the 
physical amplitudes. Although with large statis¬ 
tical and systematic errors, both groups demon¬ 
strated that a controlled numerical signal can be 
obtained for these matrix elements. The system- 
atics of these important results can be reduced 
by using fermions with an exact chiral symme¬ 
try, lighter quark masses and by considering the 
effective Hamiltonian with a dynamical charm. 

A different avenue is being followed in Ref. [|l6|. 
It is conceivable that the LECs of the weak chi¬ 
ral Lagrangian can be extracted by studying the 
weak interactions in the e-regime |l(J . A numeri¬ 
cal feasibility study of this approach on the lattice 
is under way. If feasible, it will be very interest¬ 
ing to compare the results of the LECs in the two 
regimes. 

For direct CP violation, both the ultraviolet 
and the infrared problems are more severe. An 
active charm does not mitigate the ultraviolet 
renormalization, and divergent power subtrac¬ 
tions are necessary. The cancellation between two 
large competing contributions from Q§ and Qs 


renders e'/e very sensitive to higher order correc¬ 
tions in chiral perturbation theory. Leading-order 
terms may not be sufficient to reach a reliable pre¬ 
diction in the Standard Model |93 |,[m). 

7.2. K°—K° mixing: e 
By using the OPE, the AS = 2 effective Hamil¬ 
tonian is given by 

h as =2 = &k Cl {n) dr(/i) + h.c. , (61) 

where the expression of the Wilson coefficient 
Ci(/i) is known at NLO (9f| and the correspond¬ 
ing bare four-fermion operator 

Ox = (s7 /i P-J)(s7 A1 P_J) (62) 

is multiplicatively renormalizable. The matrix el¬ 
ement that encodes the long-distance QCD con¬ 
tributions to e 

{K°\dx^)\K°) = \F 2 K M 2 K B K ^) (63) 

has been computed with overlap fermions in the 
last year (3gH)tj|. A plain overlap action has been 
used in Ref. |38f for a lattice of linear exten¬ 
sion L ~ 1.5 fm, with a spacing a ~ 0.093 frn 
and for degenerate light quark masses in the 
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range m s /2 < m < m s . The RI/MOM non- 
perturbative renormalization procedure has been 
implemented to compute the logarithmic diver¬ 
gent renormalization constant. NNC-HYP over¬ 
lap fermions have been used in Ref. |k| for a 
lattice of linear extension L — 1.5 fm, with a 
spacing a ~ 0.125 fm and for degenerate light 
quark masses in the range m s < m < 2.5 m s . 
The operator has been renormalized using one- 
loop perturbation theory. The results of the two 
computations are in very good agreement in the 
common range of simulated masses, as shown in 
Fig. | 

In Ref. Q the results have been slightly ex¬ 
trapolated to the physical point by using the func¬ 
tional form 

. * ( i-3(fD 2 .o g( M) 

+ O 4 ) <“> 

with F = FP hys , while a linear extrapolation has 
been performed in j9(J. Including the statistical 
errors only, the preliminary results 

F?f®(2GeV) = 0.61 ±0.07 @ (65) 

= 0.66 ±0.04 @ (66) 


are in very good agreement. More studies are 
needed to properly assess the magnitude of the 
various systematic errors. 

A comparison with other determinations ob¬ 
tained with different regularizations is shown in 
the second plot of Fig. |j. Even if the statistical er¬ 
rors are large, the agreement with the continuum- 
limit world averages based on staggered results 
in Ref. |)7| is very good. Results obtained with 
domain-wall fermions with a finite fifth dimen¬ 
sion are below the overlap determinations, but 


still compatible within errors [37 101 


8. Conclusions 

Studying QCD with an exact chiral symme¬ 
try at finite lattice spacing implies many theo¬ 
retical advantages: in the ultraviolet it simplifies 
the subtraction of divergences in composite op¬ 
erators, in the infrared it allows one to simulate 


massless quarks and it provides a natural defini¬ 
tion for the topological charge and the topological 
susceptibility. 

Large-scale QCD simulations with fermions 
with exact chiral symmetry are feasible with 
known algorithms and the present generation of 
computers, at least in the quenched approxima¬ 
tion. A regime of light quark masses not ac¬ 
cessible with standard fermions has already been 
reached. 

First results for the meson spectrum, light 
quark masses, the chiral condensate and Bk are 
in good agreement with previous determinations 
and suggest moderate discretization errors for 
Neuberger’s fermions. 

Important long-standing problems such as the 
A I = 1/2 rule in K —> tttt decays are greatly sim¬ 
plified and can be attacked. A combined use of 
chiral perturbation theory and fermions with ex¬ 
act chiral symmetry may lead to important new 
phenomenological informations in the next few 
years. 
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